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An  Evaluation  of  Methods  for 
Producing  Desired  Colors  on 
CRT  Monitors 


Three  experiments  were  performed  to  evaluate  methods  for 
generating  colors  having  desired  CIE  chromaticity  coor¬ 
dinates  and  luminances  on  CRT  monitors.  Methods  based 
on  piecewise  linear  interpolation  were  found  to  yield  better 
accuracy  than  others  based  on  linear  regression.  Compar¬ 
isons  between  two  monitors  revealed  no  substantial  differ¬ 
ences  in  accuracy,  but  the  monitor  s  state  of  adjustment 
was  found  to  have  important  effects.  Evidence  was  also 
obtained  concerning  the  number  of  CRT  measurement  points 
needed  to  obtain  maximum  accuracy  from  the  interpolation- 
based  methods.  The  methods'  accuracies  are  characterized 
and  avenues  for  improvement  are  suggested. 

Introduction 

An  increasing  number  of  vision  researchers  are  using  color 
CRT  monitors,  driven  by  digital  graphics  (DG)  systems,  to 
produce  color  stimuli.  Sometimes,  the  colors  must  have 
specific  CIF.  chromaticity  coordinates  and  luminances.  In 
such  cases,  these  persons  face  a  common  problem:  devel¬ 
oping  a  method  for  generating  specific  colors  on  a  DG/color 
monitor  system. 

If  only  a  few  colors  are  needed,  this  problem  can  be 
solved  via  trial -and-error.  Otherwise,  one  must  first  measure 
the  monitor’s  colorimetric  behavior  as  a  function  of  the 
values  loaded  into  the  DG  system's  digital-to-analog  con¬ 
vertors  (DACs,  which  output  the  voltages  that  modulate  the 
monitor’s  red,  green,  and  blue  electron  guns)  and  then  use 
the  resulting  data  to  develop  equations  that  predict  the  DAC 
values  needed  to  produce  specific  colors.  Many  methods 
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can  be  devised  for  accomplishing  these  two  tasks,  some  of 
which  have  been  reported'^  and  others  not.  Their  relative 
merits  are  unclear,  however,  because  little  effort  has  been 
made  to  validate  the  methods  and  quantify  their  predictive 
accuracies.  The  present  experiments  compared  the  accura¬ 
cies  of  several  methods  for  calculating  the  DAC  values  that 
will  produce  desired  chromaticity  coordinates  and  lumi¬ 
nances  on  a  DG/monilor  system. 

Seven  predictive  models,  representing  three  classes  of 
equations  (polynomials,  power  functions,  and  exponentials) 
and  two  modeling  approaches  (regression  and  linear  inter¬ 
polation),  were  compared.  All  seven  models  use  the  as¬ 
sumption  that  the  monitor’s  guns  do  not  interact.  This  has 
two  advantages.  First,  it  implies  that  the  monitor’s  colori¬ 
metric  behavior  can  be  adequately  characterized  by  mea¬ 
suring  each  gun  separately.  Each  gun  is  stepped  through  a 
series  of  DAC  values  while  the  other  two  guns  arc  set  to 
zero  and  the  monitor  is  measured  at  each  step.  The  required 
number  of  measurements  is,  thus,  much  less  than  would  be 
the  case  if  interactions  were  expected.  The  second  advantage 
is  that  the  equations  do  not  require  interaction  terms,  which 
greatly  simplifies  solving  for  DAC  values. 

Six  of  the  models  also  assume  that  the  guns’  chromaticity 
coordinates  are  invariant.  Explaining  the  consequences  of 
this  assumption  requires  discussing  tristimulus  colorimetry. 
For  any  color  produced  by  additively  mixing  three  primar¬ 
ies,  the  relationship  between  the  primaries’  luminances  and 
the  color’s  CIE  1931  XYZ  tristimulus  values  is  given  by  the 
matrix  expression 

T  =  CY,  (1) 

where  Y  is  a  3  x  1  vector  containing  the  primaries’  lu¬ 
minances,  T  is  a  3  x  I  vector  containing  the  resulting 
tristimulus  values,  and 
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i.e.,  C  is  a  3  x  3  matrix  where  xR  is  the  red  primary's  CIE 
1931  x  chromaticity  coordinate,  etc.  Equation  (1)  implies 
that 

Y  =  C~l  T,  (3) 

i.e.,  the  luminances  needed  to  produce  a  desired  color  can 
be  calculated  by  multiplying  its  tristimulus  values  by  the 
inverse  of  C.  If  the  primaries'  chromaticity  coordinates  are 
invariant,  C  will  be  also,  which  simplifies  forming  the  ma¬ 
trix.  Equations  (l)-(3)  are  also  valid  if  CIE  1976  chromat¬ 
icity  coordinates  and  tristimulus  values  are  substituted. 

The  first  six  models  described  below  use  Eq.  (3)  to  iden¬ 
tify  the  required  primary  luminances.  They  differ  only  in 
the  ways  in  which  they  then  determine  the  DAC  values  that 
will  produce  these  luminances.  Since  measurement  signal- 
to-noisc  ratio  usually  improves  as  luminance  increases,  the 
chromaticity  coordinates  used  to  form  C  in  the  present  ex¬ 
periments  were  those  obtained  for  each  gun  at  the  largest 
DAC  value  measured  during  monitor  calibration.  The  models 
are  the  following: 

PLCC  (Piecewise  Linear  interpolation  assuming  Constant 
Chromaticity  coordinates):  This  is  the  easiest  model  to  pro¬ 
gram.  Each  gun’s  luminance  is  assumed  to  change  linearly 
below,  between,  and  above  all  measured  DAC  values.  The 
adequacy  of  this  assumption  depends,  of  course,  on  making 
the  steps  between  measured  DAC  values  small.  Our  imple¬ 
mentation  of  PLCC  also  assumes  that  luminance  is  zero  at 
zero  DAC  value. 

LIN-LIN2  (linear-linear  second-order  model):  This  ap¬ 
proach  models  the  relationship  between  DAC  value  and 
luminance  for  each  gun  using  the  second-order  polynomial 
Y  -  a  +  bD  +  cD2,  where  D  is  the  gun’s  DAC  value.  Y 
is  the  resulting  luminance,  and  a.  b.  and  c  are  coefficients 
obtained  via  linear  regression.  This  model  captures  the  cur¬ 
vilinear  relationship  between  DAC  value  and  luminance.  It 
is  more  complicated  than  PLCC  because  the  software  must 
perform  second-order  linear  regression  and.  when  solving 
for  DAC  values,  must  be  able  to  decide  which  solution  is 
appropriate.  (Note  that  both  may  be  positive.) 

LOG-LOG:  The  relationship  between  DAC  voltage  and 
luminance  is  often  said  to  approximate  a  power  function, 
which  implies  that  equations  of  the  form  log  Y  =  a  +  b 
log  D  should  work.  This  also  requires  a  regression  routine, 
but  this  routine  and  the  DAC-value  solution  algorithm  arc 
simplified  by  the  use  of  a  first-order  model.  The  algorithm 
must,  however,  avoid  attempting  to  take  the  logarithm  of 
zero  or  a  negative  number.  Negative  values  for  Y  and/or  D 
can  arise  if  the  desired  color  lies  outside  the  monitor’s  color 
gamut. 

LOG-LOG2:  Cowan12  has  reported  that  LOG-LOG  does 
not  always  yield  accurate  results  and  has  recommended  a 
second-order  version  of  the  form  log  Y  =  a  +  b  log 
D  +  c(log  D)2.  This  has  the  same  disadvantages  as  LLN- 
LLN2  and  LOG-LOG. 

LOG-LIN:  One  can  fit  an  exponential  model  using  equa¬ 


tions  of  the  form  log  Y  =  a  +  bD.  This  requires  only  first- 
order  regression  and  has  only  one  root,  which  allows  easy 
solution  for  DAC  values. 

LOG-LIN2:  One  can  also  use  a  second-order  version  of 
IXXi-LIN,  i.e..  log  Y  =  a  +  bD  +  cD2,  which  has  the 
same  disadvantages  as  LOG-LOG2. 

PLVC  (Piecewise  Linear  interpolation  assuming  Variable 
Chromaticity  coordinates):  This  model  was  introduced  by 
Farley  and  Gutmann2  and  is  the  same  as  PLCC  except  the 
guns’  chromaticity  coordinates  arc  assumed  to  change  as  a 
function  of  luminance.  The  coordinates  obtained  at  each 
DAC'  value  during  monitor  calibration  are  considered  to  be 
correct.  They  are  assumed  to  change  linearly  at  intermediate 
values.  This  model  requires  interpolating  for  all  three  DAC 
values  simultaneously,  which  is  more  complicated  than  us¬ 
ing  PLCC.  Methods  for  accomplishing  this  are  discussed  in 
Appendix  A. 

Method 

Apparatus 

The  apparatus  consisted  of  a  minicomputer,  digital  image 
processor,  two  19-inch  (48-cm)  color  monitors  (Tektronix 
690-SR  and  Aydin  8830  Patriot)  with  0.31 -mm  mask  pitch, 
and  a  computer-controlled  speclroradiometer.  The  Tektron¬ 
ix  contains  a  delta-gun  CRT  and  has  achieved  an  especially 
favorable  reputation  among  some  vision  researchers.  The 
Aydin  contains  an  inline-gun  CR  T  and  is  representative  of 
many  commonly  used  units.  The  image  processor  contains 
10-bit  DACs,  i.e..  each  DAC  maps  integers  ranging  from 
0-1023  onto  its  output  voltage  range  (i.e.,  0-1  VDC). 

Basic  Procedure 

For  each  replication  of  an  experiment,  the  monitor  was 
calibrated  by  measuring  each  gun  spcctroradiometrically  at 
several  equally  spaced  DAC  values.  For  each  gun,  the  first 
DAC  value  was  the  smallest  that  produced  reliable  mea¬ 
surements  and  the  last  was  the  largest  that  could  be  loaded 
into  the  DACs  (1023).  (The  lowest  luminance  that  can  be 
reliably  measured  by  our  spectroradiometer  is  approximately 
0.02.  0.10,  and  0.04  cd/m2  for  the  red.  green,  and  blue 
guns,  respectively.  The  corresponding  DAC  values  vary 
according  to  the  monitor  and  its  adjustment.)  The  resulting 
calibration  file  was  then  used  to  calculate  the  DAC  values 
needed  to  produce  a  standard  set  of  colors,  using  various 
models.  The  standard  color  set  consisted  of  28  pairs  of 
chromaticity  coordinates  (shown  in  Figure  1)  that  fully  sam¬ 
pled  the  monitors'  chromatic  gamuts.  They  were  replicated 
at  10  and  30  cd/m2. 

Each  set  of  values  was  loaded  into  the  DACs  and  the 
monitor  was  measured  to  determine  the  resulting  chromati¬ 
city  coordinates  and  luminance.  Two  dependent  measures 
were  calculated  for  each  observation:  (1)  absolute  value  of 
the  percent  luminance  error;  and  (2)  chromatic  error,  ex¬ 
pressed  as  distance  on  the  CIE  1976  uniform  chromaticity 
scale  (CCS)  diagram.  The  order  of  the  measurements  was 
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FIG.  1  CIE  1976  UCS  diagram  showing  the  28  chromaticities  used  in  the  experiments. 


randomized  for  each  replication.  Two  monitors  were  tested, 
to  check  the  generalizability  of  the  results. 

Two  methods  of  adjusting  the  monitors  were  also  tested, 
to  evaluate  the  effects  of  monitor  setup.  The  standard  setup 
involved  adjusting  the  monitor’s  white  point  to  the  chro- 
maticity  coordinates  of  CIF.  standard  illuminant  Du  while 
maximizing  the  blue  gun’s  luminance.*  This  is  fairly  con¬ 
ventional  and  has  the  advantage  of  tending  to  maximize  the 
DG/monitor  system’s  colorimetric  gamut.  The  alternate 
("equiluminous")  setup  involved  setting  each  gun's  maxi¬ 
mum  luminance  at  approximately  35  cd/m7,  which  is  roughly 
the  highest  possible  luminance  from  our  monitors’  blue  guns. 
This  sacrifices  gamut  but  yields  a  consequent  improvement 
in  colorimetric  resolution,  because  the  range  of  possible 
DAC  values  is  mapped  onto  a  smaller  luminance  range  for 
the  red  and  green  guns,  i.e.,  the  red  and  green  DAC  values 
arc  mapped  more  densely. 

Experiment  1 

The  first  experiment  compared  the  seven  models  to  identify 
the  most  accurate  ones.  Monitor  calibration  was  performed 
as  described  earlier.  In  this  case,  wc  measured  16  DAC 


“After  completing  the  standard  setup  for  the  Aydin.  we  found  that  the 
maximum  luminance  of  the  Aydin \s  blue  gun  was  roughly  20  cd/m:.  This 
was  insufficient,  because  our  color  set  requires  luminances  as  high  a>  30 
cd/m2  from  each  gun.  Wc  could  not  get  more  luminance  from  the  red  and 
green  guns,  so  we  increased  the  blue  gun’s  gain  while  leaving  the  red  and 
green  guns  alone.  This  means  the  Aydin's  “white  point"  was  actually  a 
desaturated  blue.  The  resulting  CIE  1976  chromaticity  coordinates  were 
u'  =  0.197.  v’  =  0.401,  which  corresponds  to  a  correlated  color  tem¬ 
perature  of  34.259  K. 


values  per  gun.  Two  replications  were  performed  on  each 
monitor. 

A  full-factorial  analysis  of  variance  (ANOVA)  was 
performed  for  both  dependent  measures,  using  Method, 
Monitor,  Setup,  Luminance.  Chromaticity,  and  Replica¬ 
tion  as  the  main  effects.  Statistically  significant  effects 
(significance  was  assessed  using  alpha  =  0.05  in  all 
analyses  herein)  account  for  99.5  and  97.2%  of  the  vari¬ 
ance  in  luminance  error  and  chromatic  error,  respectively. 
These  results  reflect  the  fact  that  there  was  little  variation 
across  the  two  replications,  so  nearly  every  effect  is  sig¬ 
nificant. 

Most  of  the  effects  involving  Monitor  arc  significant,  but 
they  account  for  little  of  the  variance.  Overall,  the  Aydin 
yielded  slightly  better  accuracy  than  the  Tektronix  (average 
absolute  luminance  error  of  19.8  vs.  21.4%  and  average 
chromatic  error  of  0.014  vs.  0.015).  These  differences  have 
little  practical  importance. 

Table  I  shows  the  mean  errors  from  the  Method  x  Setup 
x  Luminance  interaction.  It  provides  a  satisfactory  sum¬ 
mary  of  the  most  important  results,  although  it  must  be 
noted  that  effects  involving  Chromaticity  account  for  much 
of  the  variation  in  chromatic  error.  Critical  differences  for 
luminance  error  and  chromatic  error  in  Table  I  were  cal¬ 
culated  using  Fisher’s  Least  Significant  Difference  (LSD) 
post  hoc  paired-comparison  procedure.  This  yielded  val¬ 
ues  of  0.8%  and  0.001,  respectively.  Using  these  criteria, 
it  appears  that  most  of  the  differences  in  Table  I  are  sta¬ 
tistically  reliable.  Therefore,  the  models  based  on  linear 
regression  are  clearly  inferior  to  PLVC  and  can  be  elim¬ 
inated.  The  same  argument  could  be  made  concerning 
PLCC,  but  it  is  the  second  most  accurate  model,  is  much 
simpler  than  PLVC.  and  its  accuracy  seems  to  be  at  least 
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TABLE  I.  Mean  errors  from  Experiment  1. 


Model 

Standard  setup 

Equiluminous  setup 

10  cd/m2 

30  edm* 

10  cd/nr* 

30  cd/m* 

l%vi 

W  V' 

l%yi 

u'  V 

l%T| 

u‘  V* 

m 

u'  V' 

PLCC 

1.9 

0  008 

2.2 

0.005 

5.6 

0.011 

4.7 

PLVC 

2.0 

0.004 

2.2 

0.003 

4.5 

0.006 

4.2 

LIN-LIN2 

10.1 

0.017 

3.7 

0.008 

9.7 

0.013 

4.1 

LOG-LOG 

46.8 

0.015 

22.3 

0.017 

49.5 

0  013 

10.6 

0.027 

LOG-LOG2 

8.2 

0.017 

7.0 

0.010 

8.9 

0.022 

11.4 

0.012 

LOG-LIN 

108.5 

0.020 

693 

0.026 

71.0 

0.015 

18.2 

LOG-LIN2 

44.6 

0.024 

13.6 

0.028 

16.4 

0.025 

15.5 

0.023 

marginally  acceptable  for  some  applications.  Therefore, 
PLCC  and  PLVC  were  retained  for  further  study. 

Table  11  shows  the  average  R2  values  obtained  from  the 
regression-based  models  for  each  monitor  setup  during  the 
experiment.  It  can  be  seen  that  their  rank  order  corresponds 
fairly  well  with  the  models’  relative  performance.  However, 
the  models'  accuracies  are  lower  than  might  be  expected, 
given  the  R2 s  Harley  and  Gutmann7  and  Neri4  have  also 
noted  that  R2 s  which  would  ordinarily  he  regarded  as  very- 
good  do  not  necessarily  assure  high  accuracy. 

Better  understanding  of  the  regression-based  models’  per¬ 
formance  can  be  obtained  by  plotting  their  regression  lines. 
Figure  2  shows  an  example,  derived  from  an  Aydin  blue 
gun  calibration.  Initially,  one  notices  that  UN-LIN2  pro¬ 
vides  a  much  better  fit  than  the  others.  In  fact,  everything 
else  is  rather  poor,  despite  the  R2 s.  This  is  because  the  others 
were  fit  (and  the  R2 s  obtained)  in  transformed  domains.  Of 
particular  interest  is  the  fact  that  the  error  which  was  min¬ 
imized  in  these  cases  was  the  log  luminance  error.  Coasid¬ 
ering  the  dependent  measures  used  in  Experiment  1,  one 
might  expect  this  to  yield  better  results  than  minimizing 
luminance  error.  However,  given  LIN-LIN2’s  accuracy  rel¬ 
ative  to  the  other  models’,  this  expectation  is  apparently 
erroneous. 

It  is  possible  to  fit  the  other  models  without  transforming 
the  data.  This,  however,  requires  nonlinear  regression  meth¬ 
ods,  which  means  one  must  either  obtain  the  coefficients 
separately  or  add  significant  complexity  to  the  DAC-valuc 
solution  routine.  To  investigate  the  promise  of  this  ap¬ 
proach,  the  Statistical  Analysis  System’s5  Nl.IN  procedure 
(multivariate  secant  method)  was  used  to  compute  coeffi¬ 
cients  for  the  data  shown  in  Fig.  2,  using  the  nontrans formed 


TABLE  II.  Mean  R2  values  from  regression  models. 


Model 

Standard 

Setup 

Equiluminous 

LIN-LIN2 

0.999 

0.999 

LOG-LOG 

0.957 

0.909 

LOG-LOG2 

0.993 

0.982 

LOG-LIN 

0.800 

0.834 

LOG-LIN2 

0.940 

0.957 
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equivalents  of  the  LOG-LOG,  LOG-LOG2,  LOG-LIN,  and 
IjOG-1  JN2  models.**  The  resulting  R2  for  LOG-LIN  is  0.998, 
but  the  visual  fit  is  not  very  good.  The  other  models  yield 
R2 s  better  than  0.999  and  provide  good  visual  fit.  Therefore, 
it  seems  justifiable  to  explore  them  further,  using  this  ap¬ 
proach.  Whether  this  will  yield  better  accuracy  than  PLCC 
(or  even  LIN-LIN2)  is  an  open  question,  but  Neri’s6  report 
dial  even  sixth-order  polynomials  arc  inferior  to  PLCC  sug¬ 
gests  otherwise. 

For  PLCC  and  PLVC,  the  standard  setup  yielded  better 
accuracy  than  the  equiluminous  setup.  Therefore,  the  re¬ 
mainder  of  this  article  only  reports  results  obtained  using 
the  standard  setup.  (Earlier  experiments  using  the  equilu¬ 
minous  setup  have  also  yielded  results  for  PLCC  and  PLVC 
which  are  inferior  to  the  standard  setup  results  reported 
here.7) 

An  attempt  was  made  to  ascertain  the  reason  for  the 
standard  setup’s  superiority,  lhe  calibration  files  were  used 
to  create  plots  of  luminance  as  a  function  of  DAC  value  for 
each  gun,  but  this  revealed  no  obvious  irregularities  for  the 
equiluminous  setup.  That  is.  there  is  no  evidence  that  re¬ 
ducing  the  red  and  green  maximum  luminances  (i.c.,  re¬ 
ducing  gain)  caused  the  guns  to  behave  peculiarly. 

The  only  noteworthy  difference  is  that  the  red  and  green 
luminances  rise  much  more  rapidly  for  the  standard  setup, 
which  means  the  smallest  measurable  red  and  green  DAC 
values  are  smaller.  For  example,  for  the  Tektronix  with  the 
equiluminous  setup,  the  smallest  measurable  red  and  green 
DAC  values  were  258  and  303,  respectively,  whereas  they 
were  both  123  with  the  standard  setup.  Since  the  maximum 
red  and  green  luminances  were  much  higher  with  the  stan¬ 
dard  setup,  though,  it  appears  that,  contrary  to  expectation, 
the  density  of  the  mapping  of  DAC  value  onto  luminance 
was  not  important.  Instead,  the  important  factor  may  be  that 
the  standard  setup  lead  to  less  extrapolation  below  the  mea¬ 
sured  range  when  solving  for  DAC  values.  This  docs  not 
mean  that  mapping  density  is  never  of  consequence.  It  may 
simply  be  that  our  10-bit  DACs  provided  more  than  was 
needed. 


-These  arc.  respectively:  Y  =  aD*.  Y  =  aD'^e'°tD\  Y  =  lO***0’. 
and  Y  - 
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FIG.  2.  Luminance  as  a  function  of  DAC  value  for  Aydin  blue  gun.  as  measured  (individual  crosses)  and  as  predicted  by 
each  linear  regression  model  (solid  and  broken  lines). 


Experiment  2 

The  next  step  was  to  perform  a  more  thorough  evaluation 
of  PLCC  and  PLVC.  This  was  accomplished  by  expanding 
the  standard  color  set  to  include  a  replication  at  3  cd/m2  and 
by  testing  the  effect  (which  will  be  termed  "resolution”)  of 
the  number  of  measurements  contained  in  the  monitor  cal¬ 
ibration  file.  The  latter  test  involved.  (1)  measuring  32  equally 
spaced  DAC  values  per  gun  during  calibration,  thereby  ob¬ 
taining  a  32-point  calibration  file;  (2)  deleting  every  other 
measurement  to  obtain  a  16-point  file;  (3)  repeating  to  obtain 
an  8-point  file;  and  (4)  using  PLCC  and  PLVC  to  calculate 
DAC  values,  using  each  file.  Two  replications  of  the  ex¬ 
periment  were  performed  on  each  monitor. 

ANOVAs  were  performed  as  before,  but  with  the  main 
effect  of  Resolution  substituted  for  Setup.  Once  again,  nearly 
evcr>-  effect  is  statistically  significant,  accounting  for  a  total 
of  98.9  and  96.9%  of  the  variance  in  luminance  error  and 
chromatic  error ,  respectively. 

The  most  important  results  arc  summarized  in  Table  III. 
which  shows  the  Luminance  x  Resolution  x  Method  in¬ 
teraction.  The  LSD  critical  differences  for  luminance  error 
and  chromatic  error  in  this  table  arc  0.8%  and  0.001.  re¬ 
spectively. 

There  are  several  points  worth  making  here.  First,  16- 
point  resolution  is  clearly  superior  to  8-point  and  seems  to 


provide  as  much  accuracy  as  can  be  obtained.  Indeed,  it  is 
even  superior  to  32-point  resolution,  at  3  cd/m2.  The  effect 
of  resolution  seems  to  diminish  as  luminance  increases, 
though  perhaps  because  computations  for  the  more  lu¬ 
minous  colors  tend  to  involve  the  upper  range  of  the  DAC 
value  versus  luminance  functions,  where  there  is  less  cur¬ 
vature. 

Second,  accuracy  improves  as  luminance  increases.  Lu¬ 
minance  error  seems  to  stabilize  somewhere  between  3  and 
10  cd/m2.  PL  VC's  chromatic  error  behaves  similarly,  but 
PLCC’s  shows  substantial  improvement  between  10  and  30 
cd/m2. 

Third.  PLCC  and  PLVC  yield  comparable  luminance  er- 


TABLE  III.  Mean  errors  from  Experiment  2. 

PLCC  PLVC 


Lum 

Res 

!%V| 

u'  V* 

l%7| 

u‘  v' 

3 

8 

20.3 

0042 

19.8 

0.044 

16 

4.2 

0  015 

4.9 

0.008 

32 

6.2 

0.018 

7.3 

0  009 

10 

8 

5.4 

0.010 

5.1 

0.007 

16 

2.0 

0.009 

1.9 

0.004 

32 

2.3 

0.009 

2.3 

0.004 

30 

8 

2.7 

0.006 

2.7 

0.004 

16 

2.4 

0.005 

2.2 

0.003 

32 

2.2 

0.005 

2.2 

0.003 
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FIG.  3.  CIE  1976  UCS  diagram  inset  showing  mean  absolute  percent  luminance  error  (upper  values)  and  mean  chromatic 
error  (lower  values)  from  Experiment  2  for  Aydin  monitor  using  16-point  calibration  resolution.  PLCC  at  3  cd/m2. 


ror,  but  PLVC  clearly  yields  less  chromatic  error,  although  though  cither  model  can  yield  good  mean  accuracy  plus 

the  differences  diminish  as  luminance  increases.  These  find-  excellent  accuracy  for  specific  colors,  they  can  also  yield 

ings  are  logical  because  (1)  PLCC  and  PLVC  model  gun  substantial  error  in  some  cases.  The  colors  for  which  error 

luminance  in  the  same  way,  (2)  shifts  in  the  guns'  measured  is  high  arc  not  always  the  same  on  the  two  monitors,  but 

chromaticity  coordinates  occur  mainly  at  low  luminances.  there  are  many  similarities. 

and  (3)  PLCC's  estimates  of  gun  chromaticity  are  derived  Effects  involving  Monitor  account  for  greater  proportions 
from  the  maximum  luminances.  Therefore.  PLVC's  DAC  of  the  variance  than  in  Experiment  I .  For  PLVC.  Monitor 

value  solutions  approach  PLCC’s,  as  luminance  increases.  makes  little  overall  difference.  For  PLCC.  however,  chro- 

FinaJly,  it  is  worth  pointing  out  that  the  values  at  10  and  matic  accuracy  tends  to  be  better  on  the  Aydin  (average  of 
30  cd/m2  for  16-point  resolution  are  very  similar  to  the  0.011  vs.  0.015).  although  the  differences  diminish  as  lu- 

comparable  values  in  Table  I.  This  provides  added  confi-  minancc  increases.  These  trends  are  also  evident  in  the 

dence  in  the  repeatability  of  our  results.  standard-setup  results  from  Experiment  1,  but  not  in  the 

The  effects  of  Chromaticity  are  illustrated  in  Figs.  3-8.  equiluminous-setup  results. 

The  figures  show  chromatic  and  luminance  error  for  PLCC  We  believe  the  effects  of  Monitor  arc  due  mainly  to 

and  PLVC  at  3,  10,  and  30  cd/m2,  averaged  across  the  two  changes  in  the  Tektronix’  chromatic  gamut  which  occurred 

replications,  for  the  Aydin  using  16-point  resolution.  They  during  the  course  of  the  experiments.  Its  CRT  was  replaced 
are  intended  to  illustrate  the  range  of  accuracies  one  might  shortly  before  we  started  data  collection,  and  the  equilu- 

expect  to  achieve.  They  also  underscore  the  fact  that,  al-  minous-setup  condition  was  run  first.  Our  chromaticity  set 
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FIG.  4.  Same  as  Fig  3.  PLCC  at  10  cd/m7. 

whs  originally  designed  to  fall  within  both  monitors*  gamuts  Experiment  3 
at  all  luminances,*  but  apparently  deviated  from  this  cri-  ... 

.erion  as  .he  Tektronix  CRT  burned  in.  The  resulting  prob-  A  ,h^1  cJxPc"men'  »«  pertbrmed  to  .es.  PLCC  and  PI.VC 
lems.  which  manifest  themselves  most  obviously  in  the  font.  al  1  °°  cdA"  Bctause  of  l,m"a,l0ns  on  cach  Sun's  max‘- 
of  increased  chn.matic  error  a.  the  gamut  edges,  are  evident  T™  ll\m,nancc-  "  was  "»•  “c  the  complete  set 

for  both  PLCC  and  PI.VC.  PL  VC.  however,  tended  to  vield  °!  28  chro™llclI'cs  >»  '«'•  'merefore.  a  subset  con- 

lower  chromatic  error  elsewhere  for  the  Tektronix  'and.  of  ehroma.ic.tcs  I.  2  3.  8.  9.  10.  14.  IS.  16.  19. 

therefore,  its  average  error  was  comparable  across  Monitor.  and  20  (s“  Fl*'  '>  was  Tw0  replications  of  the  ex- 

penment  were  performed  on  cach  monitor  using  16-point 
calibration  resolution. 


•Allhough  a  monitor's  chromatic  gamut  is  typically  portrayed  as  a  sim¬ 
ple.  invariant  triangle  on  the  chromaticity  diagram,  it  is  actually  much 
more  complex.  One  reason  is  thut  the  guns  have  luminance  limitations,  so 
the  gamut  changes  shape  and  collapses  as  these  maxima  are  exceeded. 
Furthermore,  the  guns'  measured  chromaticity  coordinates  vary  with  lu¬ 
minance.  Farley*  and  Post*  have  suggested  (hat  this  is  due  primarily  to 
signnl-to-noi&c  ratio  problems  in  the  measuring  device.  It  is  also  likely, 
though,  that  phosphor  persistence,  unintentional  excitation  by  room  lights 
and  reflections  off  the  faceplate,  imperfect  purity  adjustment,  shadowmask 
motion,  and  stray  electrons  play  a  role. 


TABLE  IV.  Mean  errors  from  Experiment  3. 


Lum 

PLCC 

PLVC 

%Y 

|%Yi 

u'  V 

%Y 

i%y| 

U'  V1 

3 

1.2 

3.8 

0.011 

0.8 

4.2 

0.006 

10 

0.0 

1.1 

0008 

-0.5 

1.1 

0.002 

30 

1.2 

1.9 

0.006 

0.6 

1.7 

0.002 

100 

-0  9 

2.3 

0.003 

-1.1 

2.1 

0.002 
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FIG.  5.  Same  as  Fig.  3.  PLCC  at  30  cd/m*. 


The  results  arc  summarized  in  Table  IV.  For  comparative  minance  error  shows  an  interesting  and  statistically  signif- 
purposcs,  means  for  the  same  subset  of  chromalicitics  at  3,  icant  tendency  to  go  from  being  positive  at  30  cd/m2  to 
10,  and  30  cd/m2  are  also  shown.  They  were  computed  negative  at  100  cd/m2.  This  suggests  subadditivitv  among 
using  the  data  from  Experiment  2.  In  addition,  the  column  the  guns,  possibly  caused  by  limitations  in  the  monitors* 
showing  absolute  percent  luminance  error  has  been  supple-  power  supplies. 

mented  with  another,  showing  percent  luminance  error.  The  Most  of  the  average  errors  shown  in  Table  IV  for  PLVC 
LSD  critical  differences  for  luminance  error,  absolute  lu-  are  well  within  customary  needs.  Indeed,  we  wondered 

minance  error,  and  chromatic  error  in  Table  IV  arc  1.3%,  whether  substantially  better  results  are  possible,  given  the 

1.0%,  and  0.001.  respectively.  inherent  variability  of  the  DG/monitor/spectroradiometer 

The  results  are  mostly  in  keeping  with  expectations  based  systems.  Therefore,  we  assessed  this  variability, 

on  the  preceding  findings.  PLVC’s  chromatic  error  stabilizes  PLVC  was  used  to  calculate  Aydin  DAC  values  for  the 
between  3  and  10  cd/m2  and  remains  stable  up  to  100  chromaticity  subset  at  all  four  luminances.  The  DAC  values 

cd/m2.  PLCC's  chromatic  error  continues  improving  up  to  were  loaded  and  measured  five  times  in  random  order  and 

100  cd/m2,  at  which  point  it  is  no  longer  reliably  worse  than  the  average  luminance  and  chromaticity  coordinates  ob- 

PLVC's.  Absolute  luminance  error  docs  not  differ  reliably  tained  for  each  set  of  DAC  values  was  calculated.  We  then 

across  models  at  any  luminance.  However,  in  both  cases,  computed  the  difference  between  each  measurement  and  the 

it  is  reliably  worse  at  100  cd/m2  than  at  10  cd/m2,  and  is  associated  average  values.  Specifically,  we  calculated  per- 

intermediale  at  30  cd/m2.  Apparently,  it  reaches  a  minimum  cent  luminance  error  and  error  in  CIE  u‘  and  v\  Next,  the 

somewhere  near  10  cd/m2  and  then  slowly  increases.  Lu-  standard  deviation  of  these  errors  was  computed  for  each 
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PI  .VC’s  chromatic  error  is  roughly  ten  times  the  standard 
deviations,  which  are  comparable  with  the  smallest  steps  in 
chromaticity  we  can  lake,  given  our  10-bit  DACs.  This 
indicates  that  substantial  reductions  in  chromatic  error  can 
probably  be  achieved  at  all  luminances. 

There  are  two  other  points  worth  making  about  Table  V. 
First,  chromatic  stability  tends  to  improve  as  luminance 
increases.  This  may  partially  explain  the  reduction  in  chro¬ 
matic  error  that  we  observed  as  luminance  increased.  Sec¬ 
ond,  with  the  possible  exception  of  luminance  stability  at 
3  cd/m7,  there  arc  no  pronounced  differences  between  mon¬ 
itors. 

Discussion 

Our  findings  have  several  practical  implications  for  persons 
who  wish  to  use  DG/monitor  systems  to  generate  colori- 
mctrically  calibrated  stimuli.  First,  PI.CC  and  PLVC  are 
clearly  better  than  the  other  models  considered  here  because 
they  are  more  accurate.  Choosing  between  them  requires 


FIG.  6.  Same  as  Fig.  3,  PLVC  at  3  cd-m?. 


TABLE  V.  Average  standard  deviations  from 


Experiment  3. 


Lum 

Aydin 

Tektronix 

%y 

u‘ 

V' 

%y 

u ' 

V' 

3 

1.4 

0.0006 

0.0006 

0.5 

0.0003 

0.0004 

10 

0.6 

0.0003 

0.0003 

0.5 

0.0004 

0.0005 

30 

0.7 

0.0003 

0.0003 

0.9 

0.0002 

0.0003 

100 

2.0 

0.0003 

0.0003 

1.9 

0.0002 

0.0002 

set  of  DAC  values.  (The  average  errors  are,  of  course,  zero.) 
Finally,  the  standard  deviations  were  averaged  to  yield  a 
summary  statistic  for  each  luminance  level.  This  procedure 
was  then  repeated  for  the  Tektronix. 

The  resulting  averages  are  shown  in  Table  V.  Evidently, 
system  variability  constitutes  an  increasing  proportion  of 
the  absolute  luminance  errors  shown  in  Table  IV,  as  lu¬ 
minance  increases.  Significant  luminance-error  reduction 
appears  possible  at  3  cd/m2,  but  the  potential  is  more  modest 
at  10  and  30  cd/m2  and  perhaps  negligible  at  100  cd/m7. 
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FIG.  7.  Same  as  Fig.  3,  PLVC  al  10  cd/m2. 


considering  ihe  application.  If  approximations  will  suffice  have  reached  the  same  conclusion,  if  they  had  performed 
or  only  high  luminances  arc  needed.  PLCC  is  belter  because  similar  analyses. 

it  is  easy  to  implement  and  fairly  accurate.  If  maximum  Our  second  major  finding  is  that,  for  either  model,  a 
predictive  accuracy  across  a  broad  range  of  luminances  is  conventional  monitor  setup  can  be  expected  to  yield  better 
required,  though,  PLVC  is  definitely  better.  accuracy  than  the  alternative  explored  here.  Possibly,  even 

Recently,  Cowan  and  Rowell10  and  Brainard1 1  have  prc-  better  results  can  be  obtained  by  deviating  from  our  standard 
sented  analyses  from  which  they  conclude  that  their  CRT  setup.  In  particular,  maximizing  each  gun’s  luminance  seems 

monitors  meet  (within  reasonable  limits)  the  assumption  that  worth  exploring  because  this  should  reduce  extrapolation 

the  chromaticity  coordinates  associated  with  each  gun  are  below  the  measured  range  when  solving  for  DAC  values, 
invariant.  They  have  used  the  term  "phosphor  constancy"  Third,  for  cither  model,  a  16-point  calibration  should 
to  refer  to  this  assumption  (although  there  is  no  special  yield  as  much  overall  accuracy  as  can  be  achieved.  Because 

reason  to  doubt  that  the  phosphors'  chromaticity  coordinates  the  upper  portions  of  the  DAC-value  versus  luminance  func- 

are  constant).  Clearly,  our  conclusions  regarding  the  merits  tions  are  nearly  linear,  though,  it  is  likely  that  fewer  mca- 

of  PLCC  versus  PLVC  imply  that  our  monitors  behave  surements  will  suffice  if  the  measurement  points  are  chosen 

differently.  This  apparent  contradiction  may,  however,  re-  properly.  (It  must  be  noted,  however,  that  monitors  exhibit 

llect  differences  in  analyses,  rather  than  differences  in  hard-  deceleration  at  the  top  end  if  the  gain  is  set  too  high.) 

ware.  We  have  actually  quantified  the  improvement  pro-  Therefore,  if  the  monitor  setup  is  relatively  permanent,  a 

duccd  by  allowing  for  variable  chromaticity  coordinates  and  reasonable  procedure  might  be  to  initially  perform  a  16- 

found  it  to  be  nontrivial.  The  aforementioned  authors  might  point  calibration,  examine  the  results,  and  reduce  the  num- 
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FIG.  8.  Same  as  Fig.  3,  PL  VC  at  30  cd’m2. 


her  of  measurements  at  the  upper  portions  accordingly  in 
subsequent  calibrations.  Whether  the  resulting  time  savings 
will  be  worth  this  trouble  depends  on  the  application. 

Fourth,  neither  model  is  very  accurate  at  luminances  of 
3  cd/m2  or  (presumably)  less.  This  may  be  inherent,  hut 
seems  more  probably  attributable  to  our  ability  to  accurately 
calibrate  the  monitors  at  low  luminances.  At  worst,  though, 
given  a  standard  setup  and  16-point  calibration,  PI. VC's 
accuracy  seems  to  stabilize  somewhere  between  3  and  10 
cd/m2  and  holds  reasonably  well  to  at  least  100  cd/m2. 
PLCC's  overall  accuracy  seems  to  improve  with  luminance 
up  to  roughly  100  cd/m2,  at  which  point  it  becomes  equiv¬ 
alent  with  PI.VC. 

Finally,  the  fact  that  two  very  different  monitors  yielded 
very  similar  results  suggests  that  our  conclusions  apply  to 
many,  if  not  most,  other  monitors.  Indeed,  it  appears  likely 
that  similar  accuracies  would  be  obtained  in  many  cases, 
given  comparable  measuring  conditions. 


Conclusions 

The  accuracies  obtained  here  are  acceptable  for  many  appli¬ 
cations.  but  not  all.  Substantial  improvements  are  apparently 
possible  (especially  at  the  lower  luminances)  and  greater  con¬ 
sistency  is  desirable.  Our  evidence  suggests  that  better  char¬ 
acterization  of  the  low  end  of  the  DAC  value  versus  luminance 
functions  would  yield  worthwhile  improvements. 

Further  error  reduction  may  be  possible  using  models  that 
account  for  interactions  among  the  monitor's  guns.  Unfor¬ 
tunately.  this  requires  much  more  complex  DAC  value  so¬ 
lution  algorithms  and  calibration  procedures.  In  any  case, 
though,  it  appears  that  alternate  models  must  account  for  shifts 
in  the  guns*  measured  chromaticity  coordinates  if  they  arc  to 
outperform  PLVC.  It  is  hoped  the  results  presented  here  will 
sen  e  as  useful  benchmarks  for  the  development  of  improved 
models. 

It  is  important  to  realize  that,  thus  far,  no  model  has  been 
shown  to  consistently  yield  high  predictive  accuracy.  There- 


182 


COLOR  research  and  application 


GREEN 
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RED  DAC  VALUE 

FIG.  9.  PLVC  solution  spaces  for  an  8-point  calibration  on  a  system  having  8-bit  DACs. 


fore,  for  applications  requiring  maximum  accuracy.  DAC 
value  calculations  should  not  be  trusted.  Instead,  they  should 
be  measured  and  corrected.  This  can  be  accomplished  by  a 
computerized  search  program  that  uses  the  calculations  as 
a  starting  point  and  then  repetitively  measures  the  resulting 
colors  using  a  computer-controlled  spcctroradiomctcr  and 
adjusts  the  DAC!  values  until  user-specified  colorimetric  cri¬ 
teria  arc  met.  An  algorithm  for  accomplishing  this  is  pre¬ 
sented  in  Appendix  B. 

Appendix  A:  PLVC 

Before  discussing  PLVC.  it  is  worthwhile  to  consider  al¬ 
ternate  approaches  that  also  allow  for  variable  chromaticity 
coordinates.  One  is  to  develop  regression  equations,  relating 
the  elements  of  C_l  to  T.  Given  T  and  these  equations.  C  1 
can  be  computed.  Eq.  (3)  can  be  solved  for  Y.  and  the  DAC 
values  needed  to  produce  Y  can  be  obtained  using  any  of 
the  first  six  methods  described  in  the  Introduction,  as  well 
as  others.  The  main  problem  with  this  approach  is  that  the 
guns'  measured  chromaticity  coordinates  do  not  change  lin¬ 
early  with  luminance.  Instead,  they  change  most  rapidly  at 
low  luminance  and  then  stabilize.  This  means  the  relation¬ 
ship  between  C"1  and  T  is  nonlinear  in  a  way  that  can  be 
difficult  to  model  accurately  in  a  regression  equation.  Fur¬ 
thermore.  the  calculations  needed  to  obtain  the  tables  of 


independent  and  dependent  variables  for  regression  arc  nu¬ 
merous  and  cumbersome. 

A  simpler  approach  is  to  fit  regression  equations  to  the 
calibration  data,  independently  predicting  each  tristimulus 
value  for  each  gun  as  a  function  of  DAC  value.  l*he  equa¬ 
tions  can  be  used  to  create  a  lookup  table,  showing  predicted 
tristimulus  values  for  every  possible  DAC  value  for  each 
gun.  This  table,  in  turn,  can  be  used  to  create  another  table, 
showing  predicted  tristimulus  values  for  every  possible  com¬ 
bination  of  DAC  values.  Solving  for  DAC  values,  then, 
becomes  a  relatively  trivial  matter  of  searching  the  second 
table  for  the  closest  tristimulus  values  and  reading  off  the 
associated  DAC  values. 

One  problem  with  this  approach  is  the  size  of  the  tables. 
If  eight-bit  DACs  are  used,  the  first  tabic  will  contain 
(3  x  3  x  2*=)  2304  values  and  the  second  will  contain 
(3  x  2563  =)  50.331.648.  For  ten-bit  DACs,  the  sizes 
become  9216  and  3,221.225,472,  respectively.  Thus,  the 
second  table  is  too  large  to  reside  in  memory  and  can  be 
too  large  for  a  single  disk,  even  if  one  is  willing  to  suffer 
the  time  penalty  for  performing  the  disk  reads.  Of  course, 
the  need  for  the  second  (and  even  the  first)  table  can  be 
eliminated  by  doing  the  calculations  "on  the  fly”  but.  in  any 
case,  the  time  required  to  obtain  a  solution  on  most  com¬ 
puters  will  be  measured  in  minutes,  if  not  hours,  and  there¬ 
fore  becomes  excessive  for  many  applications.  A  second 
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FIG.  10.  CIE  1976  UCS  diagram  showing  the  210  chromaticities  adjusted  each  night. 


problem  is  that  our  experimental  results  indicate  that  models 
based  on  linear  regression  are  not  very  accurate,  although 
the  incorporation  of  variable  chromaticity  coordinates  might 
help. 

PI  .VC  makes  an  easier  approach  possible.  To  explain  this 
approach,  let  us  simplify  matters  initially  by  assuming  that 
the  luminance  of  each  gun  is  proportional  with  its  DAC 
value.  Note  that  the  elements  of  C  arc  the  slopes  of  the 
guns*  tristimulus  values  as  a  function  of  Y.  We  can  now 
write  a  modified  version  of  Eq.  (1): 

T  -  SO,  (4) 

where  D  is  a  3  x  I  vector  containing  the  DAC  values,  T 
is  the  same  as  before,  and 


dXR/dDR 

dXc-dDc. 

dYR!dDR 

dYJdDc 

dZR/dDR 

dZcJdDc, 

dXjdD* 
dYh/dDH  , 
dZofdDfi 


(5) 


i.e.,  S  is  a  3  x  3  matrix  where  dXRtdDR  is  the  slope  of  the 
X  tristimulus  value  as  a  function  of  red  DAC  value,  etc. 
Thus,  S  is  equal  to  C,  multiplied  by  the  constants  of  pro¬ 
portionality  relating  Y  and  D.  This  leads  to 


D  S-1  T,  (6) 

which  is  a  modified  version  of  Eq.  (3).  The  most  important 
difference  is  that  Eq.  (6)  yields  DAC  values  directly. 

The  problem  in  actual  practice,  of  course,  is  that  Y  and 
O  arc  not  proportional  and,  furthermore,  the  guns'  chro¬ 
maticity  coordinates  vary  with  Y.  Within  the  context  of 
PLVC.  however,  Eqs.  (4)  and  (6)  can  be  treated  as  being 
valid  within  each  range  of  interpolation.  That  is,  we  assume 
that,  within  each  interpolation  range,  Y  and  D  are  linearly 
related  (this  assumption  is  also  used  for  PLCC)  and  the 
guns’  chromaticity  coordinates  change  linearly  with  D,  al¬ 
though  the  slopes  vary  from  one  range  to  another.  Thus, 
we  divide  the  DG/monitor  system’s  colorimetric  gamut  into 
many  small  solution  spaces,  each  having  its  associated  ma¬ 
trix  S,  and  pretend  the  system  is  linear  within  each  space. 
Figure  9  illustrates  this  subdivision  for  an  8-point  calibration 
on  a  system  having  8-bit  DACs. 

If  the  correct  space  is  known,  the  slopes  that  make  up  S 
can  be  readily  computed  from  the  calibration  file.  The  only 
remaining  complication  is  that  intercepts  must  be  subtracted 
from  T  before  solving  Eq.  6  and  intercepts  must  be  added 
to  D  afterwards  (just  as  is  the  case  when  using  PLCC).  This 
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is  fairly  easy,  though,  because  these  intercepts  arc  entries 
in  the  calibration  file,  anyway.  Thus,  in  actual  practice, 
Hqs.  4  and  6  become 

T  =  S(D-Do)  +  To  and  (7) 

D  =  S-*  (T-T0)  +  D0.  (8) 

respectively,  where  D0  is  a  3  x  1  vector  containing  the 
solution  space’s  intercepts  (i.e.,  the  smallest  DAC  values 
within  the  solution  space),  T0  is  a  3  x  l  vector  containing 
the  trislimulus  values  associated  with  D0.  and  the  other  terms 
are  as  before.  Note  that  the  slopes  of  all  three  tristimulus 
values  can  change  independently  from  one  S  to  another, 
thereby  accounting  for  changes  in  each  gun's  measured 
chromaticity  coordinates.  (If  the  coordinates  were  constant, 
the  slopes  of  X  and  Z  would  always  be  in  the  same  pro¬ 
portions  with  the  slopes  of  Y,  although  the  slopes  of  Y  could 
change  independently  from  one  S  to  another.) 

It  is  important  to  realize  that  every  solution  space  will 
yield  a  solution  to  tq.  (8).  How,  then,  can  the  correct  space 
be  recognized?  The  correct  space  is  the  one  that  yields  values 
for  D  that  lie  within  the  space.  If  one  or  more  values  lie 
outside,  they  have  been  obtained  via  extrapolation,  rather 
than  interpolation. 

The  more  difficult  problem  is  finding  the  correct  space. 
Because  interpolation  is  always  performed  between  the  two 
nearest  measured  points,  the  total  number  of  solution  spaces 
is  equal  to  the  cube  of  the  calibration  resolution.  Therefore, 
a  16-point  calibration  file,  for  example,  provides  4096  pos¬ 
sible  solution  spaces.  Thus  far.  we  have  been  unable  to 
devise  an  efficient  way  to  directly  identify  the  correct  space, 
given  T.  We  believe  it  is  possible  to  generate  a  lookup 
table,  but  this  does  not  appear  easy  and,  for  a  16-point 
calibration  file,  the  table  requites  (3  x  4  x  4096  *)49.152 
values  (each  entry  consisting  of  one  tristimulus  value  plus 
three  DAC  values).  Thus,  this  approach  docs  not  seem  very 
attractive.  We  welcome  suggestions  from  anyone  having 
further  insight  concerning  this  problem. 

Two  basic  search  techniques  are  possible.  The  simplest 
is  brute  force:  calculate  all  possible  solutions  until  the  correct 
space  is  encountered.  This  is  slow  and  inefficient.  A  better 
technique  is  to  use  the  results  from  each  trial  solution  as  an 
indication  of  which  space  should  be  tried  next.  This  usually 
locates  the  correct  space  very  quickly,  but  is  vulnerable  to 
infinite  loops.  We  use  the  latter  technique,  but  have  had  to 
include  code  to  trap  and  recover  from  loops. 

It  is  worthwhile  to  point  out  a  pitfall  that  can  arise  with 
either  search  technique.  If  the  relationship  between  DAC 
value  and  the  tristimulus  values  is  not  perfectly  monotonic 
in  the  calibration  file,  multiple  solutions  become  possible 
for  some  tristimulus  values.  That  is,  for  some  tristimulus 
values,  there  will  be  more  than  one  solution  space  that  yields 
values  for  D  that  lie  inside  the  space.  The  most  obvious 
way  in  which  this  can  occur  is  if  the  monitor's  luminance 
drifts  too  much  (given  the  calibration  resolution)  during 
calibration.  This  can  yield  a  measured  luminance  at  one 
DAC  value  that  is  higher  than  at  the  next  larger  measured 
DAC  value,  thereby  yielding  a  negative  slope  for  Y  (and. 


probably.  X  and  Z)  within  that  interpolation  range.  A  less 
obvious  mechanism  involves  shifts  in  the  measured  chro¬ 
maticity  coordinates.  This  can  cause  the  slope  of  X  and/or 
Z  to  become  negative,  even  though  the  slope  of  Y  may  be 
positive.  Therefore,  if  one  chooses  to  check  the  calibration 
file  for  non-monotonicity,  all  three  tristimulus  values  must 
be  checked — checking  Y  alone  is  not  sufficient. 

One  final  pitfall  is  worth  discussing.  Occasionally,  it  will 
be  found  that  no  solution  space  meets  the  criterion  for  cor¬ 
rectness.  even  though  the  color  lies  within  the  system’s 
gamut.  This  occurs  sometimes  for  colors  that  lie  along  the 
edges  separating  the  solution  spaces  and  is  due  to  rounding 
and/or  truncation  in  the  calculations.  The  basic  problem, 
then,  is  that  one  solution  space  indicates  that  the  color  lies 
just  within  an  adjacent  space,  hut  that  adjacent  space  in¬ 
dicates  that  the  color  lies  just  within  the  other.  This  is  a 
special  case  of  the  infinite  loop  problem  which  was  men¬ 
tioned  earlier  and  is  easy  to  trap.  Choosing  between  the 
solutions  Is  rather  arbitrary,  because  the  associated  DAC 
values  rarely  differ  by  more  than  I .  A  reasonable  selection 
rule  is  to  choose  the  solution  yielding  the  smaller  predicted 
colorimetric  error. 

Appendix  B:  A  “Measure  and  Adjust"  Algorithm  for 
Producing  Desired  Colors  on  Displays 

We  have  devised  a  simple  algorithm  for  the  computerized 
search  mentioned  at  the  end  of  our  article.  Our  experience 
has  been  that  it  is  efficient  and  robust: 

1 .  Calculate  the  DAC  values  needed  to  produce  the  de¬ 
sired  color.  Call  these  the  "calculated"  values. 

2.  Ix>ad  the  calculated  values,  measure  the  display,  and 
check  the  colorimetric  tolerance.  If  the  color  is  within 
tolerance,  exit  with  the  calculated  values.  Otherwise, 
proceed  to  Step  3. 

3.  Calculate  the  DAC  values  needed  to  produce  the  color 
that  appeared  on  the  display,  using  the  same  com¬ 
putational  method  as  in  Step  I .  Call  these  the  “mea¬ 
sured"  values. 

4.  Subtract  the  measured  values  from  the  calculated  val¬ 
ues.  Then,  add  the  differences  to  the  loaded  values, 
yielding  the  “corrected"  values.  In  other  words: 

Corrected  =  Loaded  +  (Calculated  -  Measured) 

5.  Load  the  corrected  values,  measure  the  display,  and 
check  the  tolerance.  If  the  color  is  within  tolerance, 
exit  with  the  corrected  values.  Otherwise,  return  to 
Step  3. 

F.xamplc:  If  the  calculated  DAC  values  are  100,  100,  UK), 
but  the  color  they  produce  on  the  display  is  out  of  tolerance 
and  is  what  would  have  been  expected  for  1 10,  110,  110, 
the  corrected  values  become  (100  +  100  -  1 10  =)  90, 
90,  90.  If  the  color  produced  by  90.  90,  90  is  also  out  of 
tolerance  and  is  what  would  have  been  expected  for  98,  98, 
98,  the  corrected  values  become  (90  +  100  -  98  =)  92, 
92,  92.  And  so  on. 
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It  can  be  seen  that  the  algorithm  models  the  DAC-vallie 
error  as  a  simple  linear  offset  in  three  dimensions,  i.e.,  as 
though  the  primaries’  DC  offsets  (“black  levels”)  have  drifted. 
This  is,  of  course,  an  oversimplification  because  of  the 
typically  nonlinear  relationship  between  DAC  value  and 
display  luminance.  However,  if  the  original  solution  is  close 
to  being  correct,  the  linear  approximation  can  be  expected 
to  at  least  reduce  the  error,  which  means  the  approximation 
will  be  even  more  adequate  on  the  next  iteration. 

In  practice,  we  have  found  that  the  algorithm  works  well 
(at  least,  for  CRT  monitors),  even  when  the  original  solution 
errs  considerably.  For  several  years,  we  have  been  per¬ 
forming  color  vision  experiments  involving  210  colors  (see 
Fig.  10)  having  luminances  ranging  from  18  to  36  cd/m2. 
All  210  colors  arc  adjusted  each  night  to  tolerances  of 
±2.5%  in  luminance  and  0.0025  uV,  using  our  algorithm 
plus  PLVC.  On  the  first  night  following  calibration,  we 
typically  find  that:  (1)  roughly  75%  of  the  colors  are  out  of 
tolerance;  and  (2)  the  program  requires  an  average  of  slightly 
less  than  three  measurements  to  bring  each  of  these  colors 
into  tolerance.  Thereafter,  the  program  uses  the  previous 
night's  corrected  DAC  values  as  the  starting  point  and, 
again,  typically  requires  an  average  of  slightly  less  than 
three  measurements  to  bring  “bad"  colors  into  tolerance. 
This  is  near  ideal  performance  because  all  mcasurc-and- 
adjust  algorithms  require  a  minimum  of  two  measurements 
if  the  original  color  does  not  meet  tolerance. 

There  are  several  additional  points  worth  making  about 
our  algorithm.  First,  it  is  independent  of  the  DAC-valuc 
compulation  method,  as  well  as  the  method  for  assessing 
colorimetric  tolerance.  Second,  it  should  work  for  any  ex¬ 
isting  color  display  technology — not  just  CRTs.  Third, 
damping  should  be  included  in  Step  4  to  guard  against 
possible  oscillation  should  the  original  DAC  values  be  grossly 
wrong.  (We  have  found  that  a  good  maximum  permitted 
step  size  is  2‘v_4,  where  A  is  the  number  of  bits  in  each 
DAC.)  Finally,  to  avoid  infinite  loops,  two  additional  ter¬ 
mination  criteria  are  advisable: 

1.  If  the  corrected  DAC  values  lor  a  given  iteration  are 
identical  to  those  that  were  measured  in  either  the 
current  or  a  previous  iteration.  (This  indicates  that  the 
tolerances  are  loo  stringent  for  that  color,  given  the 
display's  gamut,  the  IXi/display/spcctroradiomctcr 
system's  colorimetric  resolution,  and/or  the  system's 
temporal  stability.) 


2.  If  the  number  of  measurements  exceeds  some  (gen¬ 
erous)  user-specified  limit.  (This  provides  final  pro¬ 
tection  against  gross  system  instability.) 

In  both  the  above  cases,  the  routine  should  exit  with  the  mea¬ 
sured  DAC  values  that  came  closest  to  meeting  tolerance. 
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